Unusual liquid state of hard-core bosons on the pyrochlore lattice 
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We study the physics of hard-core bosons with unfrustrated hopping (t) and nearest neighbour 
repulsion (V) on the three dimensional pyrochlore lattice. At half-filling, we demonstrate that 
the small V/t superfluid state eventually becomes unstable at large enough V/t to an unusual 
insulating state which displays no broken lattice translation symmetry. Equal time and static 
density correlators in this insulator are well-described by a mapping to electric field correlators 
in the Coulomb phase of a (7(1) lattice gauge theory, allowing us to identify this insulator with 
a U{1) fractionalized Mott insulating state. The possibility of observing this phase in suitably 
designed atom-trap experiments with ultra-cold atoms is also discussed, as are specific experimental 
signatures. 

PACS numbers: 75.10.Jm 05.30.Jp 71.27.-f a 



Much of our current understanding of the low temper- 
ature behaviour of condensed matter systems is based 
on highly successful theoretical paradigms such as Lan- 
dau's Fermi liquid theory of normal metals, Bogoliubov 
theory for superfluids, BCS theory of superconductivity 
and spin-wave theory for ferromagnets and antiferromag- 
nets However, some systems exhibit behaviour that 
falls outside of any of these standard paradigms — one 
example of this is the unconventional normal state of 
underdoped high-Tc superconductors [1, Q , while other 
examples include the cooperative paramagnetic state of 
frustrated magnets [J], and the unusual phenomenology 
of heavy fermion compounds [H]- For instance, in the un- 
derdoped normal state of high-Tc superconductors, some 
of the experimental evidence is suggestive of the fact that 
the elementary quasiparticles excitations are not spin- 1/2 
charge-e holes, but spiniess charge carriers propagating 
separately from chargeless spin carriers 0, |J. 

This has motivated much of the recent effort aimed 
at providing theoretically consistent descriptions of low 
temperature phases of matter that would display such 
spin-charge separation, or more generally, quasiparticle 
fractionalization. These developments [6] allow one to 
conclude that such exotic behaviour is indeed possible, 
and go on to provide a description of quasiparticle frac- 
tionalization in terms of an effective field theory with 
gauge symmetry p, |8j]. In this approach, fractional- 
ized quasiparticles emerge as the true low energy exci- 
tations in deconfined phases of a gauge theory (in which 
the emergent gauge force is not strong enough to bind 
the fractionalized quasiparticles into more conventional 
quanta), and can be accompanied by additional gauge 
excitations that carry energy but no spin or charge (such 
as the vortex excitation of a Z2 gauge theory Q). 

A closely related strand of activity has focused on 
the analysis of particular microscopic models in order 
to understand whether they exhibit such exotic phases 



for specific values of input parameters. This has led, for 
instance, to the construction of several different mod- 
els fill . [12I 1 which exhibit so called Z2 deconfined 
phases (the nomenclature refers to the effective gauge 
theory that affords the most 'natural' description of the 
low-energy physics) . 

One may now ask: Is there an experimental system 
which would display one of these fractionalized phases 
for a definite range of control parameters? A promising 
avenue in this regard is the physics of ultra-cold atoms in 
optical lattice potentials. Recent work has demonstrated 
that a wide variety of phenomena of interest to condensed 
matter physics can be studied by appropriately engineer- 
ing systems of ultra-cold atoms in optical potentials. For 
instance, it has been possible to provide a cold-atom real- 
ization of the superfluid-insulator transition in a bosonic 
hubbard model with on-site interactions on a cubic lat- 



tice [13|, |lj| . This has been followed by several interest- 
ing proposals for realizing fermionic and bosonic models 
with a variety of tunable interactions in different optical 
lattice geometries 15, 1^. 

In this work, we use sophisticated Quantum Monte 
Carlo (QMC) methods to provide the first confirmation 
of the existence of a U{1) fractionalized insulating phase 
that may be realized in cold-atom systems modeled by 
the Hamiltonian: 



H = J2^V{n, - l/2){n, - 1/2) - t{b\h, + hh])] 



+ Y.[U{n,- 1/2)^ - fin,] . 



(1) 



Here, rii is the particle number at sites i of a three dimen- 
sional pyrochlore lattice (Fig[T]a), bl is the corresponding 
boson creation operator, U is the on-site repulsion and V 
the nearest neighbour repulsion between bosons hopping 
(with amplitude t) on the nearest neighbour links (ij) . 
Although the pyrochlore lattice geometry we consider 



2 




FIG. 1: (color online), (a) Pyrochlore lattice and the underly- 
ing diamond lattice, (b) Ring exchange process on plaquettes 
of the diamond lattice. 



is technically challenging to realize, recent work that ap- 
peared as our study was underway provides a viable pre- 
scription for experimentally realizing such an optical lat- 
tice 



17|. Furthermore, the simplicity of the interactions 



means that they can be realized in state-of-the art cold 
atom experiments for a wide range of values of param- 
eters including the 'hard-core' limit (n^ — 0, 1) of 
very large U . We therefore focus on this hard-core limit 
in some detail here, setting t = \ and /i = in what fol- 
lows. [In this hard core limit, Eqn[l]may also be written 
in spin S' = 1/2 language via the mapping SI = ni — \l2, 

In this hard-core limit, with /i = to enforce density 
1/2 per site, the physics at small V is readily tractable: 
As the hopping t is unfrustrated, there is a stable super- 
fluid phase at small V — indeed a reasonable variational 
wavefunction for the ground state in this regime may be 
easily written down in spin language as j^*) — Y[i \Sf — 
+1/2) j. What is the low temperature state in the op- 
posite, large V limit? To answer this, we use the well- 
documented ^18] stochastic series expansion (SSE) QMC 
method (at large values of V, modifications developed 
recently !19| are crucial to maintain ergodicity — for a re- 
view, see Ref. [2(J]). 

Numerics: Most of our results are on L x L x L 
(L, the number of up pointing tetrahedra that fit in one 
side-length) samples with periodic boundary conditions 
and even L ranging from _L = 6toL = 12, and inverse 
temperature P ranging from 6 to 120 (with the largest 
(3 employed for the largest size). We use standard SSE 
estimators (isj to calculate the specific heat, the super- 
fluid stiffness ps, the bond (kinetic) energy correlations, 

and the equal time C"" (q, r = 0) = {na{q)n^' {—q)) 
and static correlators 5*"" (q, w„ = 0) = drC""^ (q, t) 
of the density nf (here a, a refer to different basis sites 
in a unit cell, and all site types (Fig. [T] (a)) are assigned 
coordinates of site- type 0). 

As is clear from Fig. [2] a), we see a distinct transi- 
tion from a superfluid state at small V, to an insulating 
state at large V for a sequence of low temperatures. This 
transition is first-order at non-zero temperature (Fig. [2] 
a), and while the first order nature is less prominent in 
lower temperature scans, a scaling analysis suggests that 
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FIG. 2: (color online), (a) Superfluid density at (3 — 30 — the 
break around V = 19.2 indicates observed hysteresis near the 
(weakly) first order transition, (b) Schematic phase diagram: 
black dots with error bars denote observed transitions, and 
red dot denotes location at which insulating phase data is 
displayed in Fig |3] and Fig|4] 
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FIG. 3: (color online). Static S""{q^,uin ~ 0) and equal time 
(^""(q, r — 0) correlators of the density in the deconfined 
phase for /3 = 30, V = 19.4. The lines are a fit to the predic- 
tions of non-compact (7(1) gauge theory on the diamond lat- 
tice as discussed in text, with fit parameters displayed above. 



the transition remains first order even in the zero temper- 
ature limit [21]. We estimate that this zero temperature 
transition is at iV/t)^ ~ 19.3 (Fig[l]b). 

In the insulator, we see absolutely no Bragg peaks that 
would correspond to spatial ordering in either the local 
density or the local bond energy. The insulator is thus, 
in this specific sense, a liquid state of matter; this is illus- 
trated in Fig [3] with several scans of density correlators 
in q space at a representative point at very low temper- 
ature above the insulating ground state. This absence of 
spatial ordering in the insulating state of an interacting 
boson system at 1/2 filling is one of our striking results, 
for such featureless insulating states are more typical of 
insulators with integer density per site. 

Interpretation: Theoretical interpretation of this strik- 
ing result is facilitated by noting that our Hamiltonian 
in this hard-core limit is closely related to that studied 
in Ref. [22: Hermele et al. considered the S* = 1/2 XXZ 
antiferromagnet on the pyrochlore lattice. By an analysis 
of a related effective model of planar rotors (with addi- 
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FIG. 4: (color online), a) The same values of sjnrj and 
^ \f^Tl fit data at different temperatures for fixed Vjt. b) 



^""(q.t^^ = 0) and /3C""(q,r = 0) for V 
tially equal at /9 = 1, but not at /3 = 30. 



19.4 are essen- 



tional terms added by hand to ensure better theoretical 
control), they argued that a deconfined phase was 
a theoretically consistent possibility in the limit of ex- 
tremely anisotropic exchange Jz 3> J_l > — however, 
since the positive sign of J_l introduces a sign problem 
in quantum Monte-Carlo treatment of such models, their 
work stopped short of making definitive statements about 
the actual phase diagram of the S' = 1/2 model. [For a 
different ejfective model of rotors that displays a U{1) 
deconfined phase, see Ref [23| 

Although our situation differs from Ref. crucially in 
the opposite sign of the hopping term, this change of sign 
does not affect 22| the arguments that make plausible 
the existence of a deconfined phase at large V/t: As the 
classical ground state for t = is macroscopically degen- 
erate (with all configurations with two particles occupy- 
ing each tetrahedron having minimum energy), the fate 
of the system at large but finite V/t is then determined 
by the structure of the effective Hamiltonian obtained to 
leading order in degenerate perturbation theory in t/V. 

As the t = ground states can be represented in terms 
of dimers living on the links of the dual diamond lat- 
tice (Fig[Ua)) subject to the constraint that two dimers 
touch each diamond lattice site, such a perturbative anal- 
ysis yields (at leading 0{t^/V'^) order) a quantum dimer 
model 22] with a 'ring-exchange' term which causes a 
'flippable' hexagon (see Fig. [1] (b)) to resonate between 
its two allowed configurations. The sign of this term is 
a matter of convention [2^ as it can be changed by an 
appropriate canonical transformation. However its struc- 
ture, and the structure of the constraint that defines the 
low energy manifold is highly reminiscent of a (compact) 
U{1) lattice gauge on the diamond latticejS^I. As the 
compact U{1) theory in three spatial dimensions admits 
a 'Coulomb' phase that mimics ordinary electrodynam- 
ics, we conclude, following Hermele et ai, that this is a 
consistent possibility at large V in our boson hubbard 
model. 

Fits: To proceed further, we note that in such a de- 
confined phase, the low energy properties are expected to 
be described by the lattice version of standard Maxwell 
electrodynamics with Hamiltonian Ti. = ^ J2<rr'> ^rr' + 



f J^oi^rr' ^ «rr')^j where the lattice curl of the vector 
potential defined on the hexagonal plaquettes of 

the dual diamond lattice (Fig[T{a)), the microscopic den- 
sity operator is related to the electric field by n^r' = 
with e a non-universal scale factor, and 7 and k are the 
emergent energy scales of this low energy description. 
[Although 7 ~ t/ — > cxD, K ~ t'^/V^ ^ in the formal 
V —^ 00, U 00 limit, their actual, renormalized values 
can be substantially different from these bare estimates.] 
To explore the implications of this ansatz for the den- 
sity correlators, it is useful to work with the correspond- 
ing imaginary time action 
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Here, f = \f>v^T is the dimensionless imaginary time vari- 
able obtained by scaling r by the typical photon energy 
Y^K7 of this artificial electrodynamics, a^r' = v^^'^a„i 
the rescaled vector potential, and a,- — v^^'^ar / ^Jlvy 
the dimensionless scalar potential (v = k/^). As den- 
sity correlators are obtained by calculating correspond- 
ing correlators of ew^/'*(9f a^r' — A^r'ar) using this ac- 
tion, it is immediately clear that this electrodynamic 
ansatz predicts C(r — 0, q) — e'^y^feq{(3y^^,q) and 

S{uj = 0,q) = e^^^=/st(q). In order to test this ansatz, 

we have calculated the functions feq{Py/Kj, q) and /st(q) 
and performed detailed fits of our data for the density 
correlators C(t = 0,q) and 5(tJ ~ 0,q). 

Our fitting procedure is quite straightforward: We first 
determine the best fit value of the scale factor Cst by 
which the function fst needs to be scaled to fit the static 
correlators S. Next, we note that the shape of feq (as a 
function of q) depends significantly on the value of the 
typical photon energy ^/kt/ that enters its first argument, 
and determine its best fit value such that feq ((S^/KFj, q) 
best mimics the shape of the corresponding equal time 
correlators Ceg(q). Finally, we determine the best fit 
value of the corresponding equal time scale factor Ceq by 
which the function f^q needs to be scaled to fit the overall 
magnitude of the equal time correlators C^q- 

Clearly, this is a very over-determined fit, since the 
same set of parameters have to fit scans of the correlators 
in the entire brillouin zone, in addition to fitting data 
at different temperatures (at fixed V/t). In addition, 
this procedure has an in-built consistency check, since 
the value of photon energy scale yTRFj can be re-obtained 
from the scale factors by noting that Ceq/cst = \f^- 

In Fig. [31 we show the results of such a fit of the static 
(zero frequency) and equal time density correlators along 
several scans in the Brillouin zone for a representative 
low temperature point at which = (similar fits work 
equally well at other low temperature points in the insu- 
lating phase) . Clearly the data fits the predictions of non- 
compact electrodynamics extremely well, with the best 
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fit values of the photon energy scale yj^rj and f^^J k,/^ 
shown in Fig. [3] (the quoted uncertainty in the best fit 
value of yTRrj also takes into account the accuracy with 
which the self-consistency condition is satisfied). Fur- 
thermore, for fixed V/t, the same parameters do indeed 
continue to fit the data as the temperature is varied (Figd] 
a)). 

Discussion: These fits are extremely convincing evi- 
dence that we have accessed the low temperature regime 
just above a Coulomb liquid ground state. Nevertheless, 
it is instructive to play devil's advocate and ask if the 
measured correlators can satisfy the predictions of non- 
compact electrodynamics to this level of accuracy if the 
system does not have a deconfined coulomb liquid ground 
state? Perhaps surprisingly, the answer is yes, but only if 
the data has been taken at temperatures such that ther- 
mal fluctuations overwhelm all quantum dynamics. 

Assume for instance that the insulator is a more con- 
ventional lattice-symmetry breaking crystal that one may 
have expected at fractional filling. If one happens to 
be at a temperature above the melting temperature of 
this putative crystal, then thermal fluctuations would 
completely overwhelm quantum effects, and the physics 
would be essentially classical. As long as the tempera- 
ture remains much smaller than V , this classical physics 
is correctly described by the classical dimer model on the 
diamond lattice, regardless of the quantum ground state. 

Now, static and equal time correlators of any quantum 
system are proportional to each other (with proportion- 
ality constant (3) in any such effectively classical regime. 
Thus, we expect dimer correlations of the classical dimer 
model to correctly describe the functional form of both 
the static and equal time correlators of the system in 
this regime. Furthermore, these classical dimer correla- 
tors are known 2j| to have precisely the same functional 
form as the static correlators of Ti. 

Regardless of the quantum ground state, we thus ex- 
pect our data for static and equal time correlation func- 
tions to be necessarily ■proportional to each other and 
match predictions of quantum electrodynamics in this 
classical regime. Is this 'trivial' mechanism responsible 
for the extremely good flts shown in FiglS]? The answer 
is clearly no: If this were the case, the static and equal 
time correlators, being proportional to each other, would 
have the same shape (as a function of q). This is clearly 
not the case for the low temperature data shown in Fig [31 
as is underscored by a comparison to data at much higher 
temperatures (Fig|4]b)), where this commonality of shape 
does become clearly visible. 

The weight of all this evidence thus allows us to con- 
clude that we are indeed seeing a Coulomb liquid state 
of matter in our simulations. What would be the best 
way to 'look' for this state of matter in a putative cold- 
atom experiment? At the most gross level, this phase is 
an incompressible insulator, with a gap to charged ex- 
citations. The distinctive difference from ordinary Mott 



insulating phases (such as those seen in the experiments 
of Ref . il3j) is the presence of a gapless neutral collective 
mode, namely the artificial photon of the U{\) gauge the- 
ory mentioned above. As we have demonstrated above, 
this neutral mode leads to characteristic dipolar struc- 
ture in the low temperature equal time and static density 
correlators. These correlations can be measured in atom- 
trap experiments by noise correlation [25| measurements 
that probe equal time correlators, and Bragg scattering 
experiments [2 a ] that probe static correlators. 
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